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An investigation on the properties of electronic states of a tight-binding Hamiltonian on the 
Apollonian network is presented. This structure, which is defined based on the Apollonian packing 
problem, has been explored both as a complex network, and as a substrate, on the top of which 
physical models can defined. The Schrodinger equation of the model, which includes only nearest 
neighbor interactions, is written in a matrix formulation. In the uniform case, the resulting Hamilto- 
nian is proportional to the adjacency matrix of the Apollonian network. The characterization of the 
electronic eigenstates is based on the properties of the spectrum, which is characterized by a very 
large degeneracy. The 2n/3 rotation symmetry of the network and large number of equivalent sites 
are reflected in all eigenstates, which are classified according to their parity. Extended and localized 
states are identified by evaluating the participation rate. Results for other two non-uniform models 
on the Apollonian network are also presented. In one case, interaction is considered to be dependent 
of the node degree, while in the other one, random on-site energies are considered. 
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I. INTRODUCTION 

Modeling physical systems on substrates without 
translation invariance provides useful insights for the un- 
derstanding of disordered systems Quite recently, 
the rapid development of complex network theory Q has 
motivated the investigation of physical models on such 
substrates, for the purpose of identifying new aspects 
that may stem from novel geometrical constructs [H, 0] ■ 
Indeed, the characterization of their physical properties 
has shown that new thermodynamic, magnetic and trans- 
port features do emerge. In several situations they 
can become quite distinct from those obtained by other 
disordered models. 

The Apollonian network |6j(AN) appears in this sce- 
nario as an interesting structure, in the sense that it 
shares several properties of complex networks but, on 
the other hand, it is defined by strict geometrical rules. 
This particular feature makes it possible to explore scale 
invariance properties together with numerical methods in 
the investigation of distinct models. Its specific geomet- 
ric features induce the emergence of non typical behavior 
in thermodynamic and magnetic properties of spin sys- 
tems @, [1] , avalanche distribution in sand pile models @ , 
and so on. 

The recent creation of a synthetic nanometer-scale 
Sierpinski hexagonal gasket, a self-similar fractal macro- 
molecule has highly motivated the area of scale in- 
variant networks. This perspective raises the question 
about electronic systems on these networks. 

The purpose of this work is to investigate tight-binding 



electronic models on the AN, focusing on the localized vs. 
extended nature of their wave functions. This characteri- 
zation follows a previous analysis of the properties of the 
eigenvalue spectrum of the adjacency matrix (AM) of the 
AN [i"l| . The AM spectrum is a very important network 
signature and, as such, its properties have been detailed 
explored within graph theory. In particular, great at- 
tention has been devoted to understanding iso-spectral 
graphs that are not isomorphic [12, EH • 



As we will show, the Hamiltonian of the uniform tight- 
binding model can be written in terms of the AM net- 
work, so that the energy eigenvalues are proportional to 
those of the AM. The previous identification of the struc- 
ture of the spectrum detailed the presence of several dis- 
tinct classes and highly degenerated subsets. The charac- 
terization of the localization properties of the wave func- 
tion, based on the evaluation of the participation rate of 
each state, follows the classification scheme of different 
sets in the AM spectrum. 



The paper is organized as follows: In Section II, we 
discuss the main features of the AN. We also introduce 
the tight-binding model and its matrix. We also indi- 
cate how the participation rate has been used to quantify 
the localization properties of the eigenstates. Section III 
brings a brief review of the most important properties of 
the uniform eigenvalue spectrum. Results for the three 
different models are presented in Section IV. Finally, Sec- 
tion V closes the paper with a summary of the our main 
important results. 
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FIG. 1: Generations n = 0, 1, and 2 of AN. The adopted site 
numbering is the one introduced in [3] 



II. TIGHT-BINDING MODELS ON AN 

AN's @, follow from the old problem of optimally 
filling a compact domain in the Euclidean plane by circles 
[TH . The solution to this problem requires the successive 
addition of maximal tangent circle into the empty re- 
gions bounded by three previously placed circles. At a 
given generation n, the network AN n is constructed by 
establishing a connection between the centers of all pairs 
of tangent circles. The solution to the packing problem 
requires the determination of the center and radius of 
each placed circle. For the AN, however, these details 
are not required, and we just need to assign the pairs 
of nodes that are connected. In this work, we consider 
the most simple situation, which arises by considering, at 
the zeroth generation n = 0, three tangent circles with 
the same radius, the centers of which occupy the cor- 
ners of a equilateral triangle (see Figure 1). The num- 
ber of sites N(n) on the network increases according to 
N(n) = (3™ + 5)/2, while the number of edges B{n) de- 
pends on n as B{n) = (3™ +1 +3)/2. In the limit n — > oo, 
B(ri)/N(n) — > 3, so that each node has, in the average, 
6 neighbors. 

The resulting AN shares several properties with other 
classes of complex network, including being scale free (the 
distribution P(k) of node degree k is a power law), small 
world (the mean minimal path (£} ~ logN, and hierar- 
chical (the clustering coefficient of individual nodes c(k) 
has a power law dependence on k). Nevertheless, it has 
own features [l6[ that precludes being put into the sets 
of networks generated by known algorithms, like those 



17] and Barabasi and 



proposed by Watts and Strogatz, 
Albert 

The model considers that an electron is primarily 
bound to a network node i with on-site energy cj, but it 
is also allowed to hop from i to a neighboring site j one 
with a probability Vi t j — Vj t i, resulting from the overlap 



integral involving the two wave functions and interaction 
potential. For a given n, the corresponding Hamiltonian 
is written as 



n)\i)eAi\ + 



E 

ihj) 



(1) 



where \i) represents a Wannier function for an electron 
at position i, and denotes pairs of sites that are first 
neighbors. 

The Hamiltonian |T]) can be written in a matrix form, 
where the diagonal elements Hi^ — e^, while off-diagonal 
elements Hij are set to Vij or 0, depending on whether 
sites i and j are neighbors or not. The AN sites are 
numbered according to a previously introduced scheme, 
which is reproduced in Figure 1 for n — 0,1, and 2. When 
we go from generation n to n + 1 , the whole n — th lat- 
tice is first squeezed into the region limited by two outer 
vertices and the site at the geometrical center of the tri- 
angle. After this first step, two identical copies of the 
squeezed lattice, rotated by 27r/3 and 47r/3, are added to 
the regions that have been emptied. In this numbering, 
two of the outer vertices always take the numbers 1 and 
2, while the third one assumes the number N(n) for gen- 
eration n. As a consequence, the central node receives 
always the number N(n — 1). The matrix representation 
of H 2 is: 
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The simple uniform model corresponds to setting = 
eo, Vi, and Vij — Vo, V(i,j). If we let E' be the eigenval- 
ues of H, we can define new energies E = [E' — £o)/Voj 
which can also be directly obtained by rewriting @ , set- 
ting eo = 0, and Vq = 1. In terms of these new vari- 
ables, the matrix ([2]) coincides exactly with the AM of 
the AN. Thus, the cigcn-cncrgics of the homogeneous 
tight-binding model can be directly obtained from the 
AM spectrum discussed previously. For the sake of a 
simple notation, we will always refer, in the forthcoming 
discussion, to the redefined values of E. Note, however, 
that even for the uniform model, the choice of the on-site 
energy eo can affect the eigenstate localization properties 
With the exception of the random on-site energy model, 
we restrict ourselves to investigate the situation in which 
£o = 0. 

Non uniform models can be defined by choosing other 
values either for the self-energies or the hopping integrals. 
In this work we consider a family of models defined by 
d = 0V«, and V(i,j) = l/z 13 . Here, z = (k l + k j )/2 is the 
average value of the node degrees k, and the value of j3 
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selects a particular element of the set. The uniform model 
corresponds to setting (3 = 0. Finally, the third model 
aims to investigate the existence of Anderson transition 
(l9| in AN is defined by setting Vri, j) = 1 for neighboring 
sites, while the t{ are chosen according to a probability 
distribution function. 

The solution to the Schrodinger equation H\^>e) = 
E\^>e) can be written in the basis of Wannier states, i.e., 
\^e) = J2i=i a E,i\i)- The site coefficients as,i are the 
components of the normalized eigenvector corresponding 
to the eigenvalue E, which are obtained from the numer- 
ical diagonalization of the Hamiltonian matrices. The 
knowledge of the vector components at and the possibil- 
ity of evaluating it as a function of the the subsequent 
generations n, suggest the use of the participation ratio £ 
to characterize the degree of localization of the electronic 
states. For normalized states, this often used measure is 
defined by 

Ei=l \°>E,i\ 

It can assume values in the range [l,iV], which corre- 
spond, respectively, to the extreme situations of com- 
pletely localized (ciea = 1 for one single value of i), and 
extended states (a^ , = 1/%/iV, Vi). Thus, it provides 
both an estimate for the localization length A, as indica- 
tions for extended states, for which £ should scale with 
N(n). Since the model is defined in the limit n — > oo, the 
localized/extended character of the Apollonian states of 
the infinite system can be inferred from the behavior of 
^E,n as function of n. 

III. SPECTRUM PROPERTIES 

The evaluation of the eigenvalues must be performed 
with the help of numeric algorithms. The spectrum of the 
uniform model has been characterized in a detailed way. 
In this particular case, the presence of many eigenvalues 
with high degree of degeneracy, and the recurrent pres- 
ence of the same eigenvalues in the successive increasing 
values of n, makes it possible to present a classification 
scheme for the eigenvalues E. According to it, for any 
value of n > 2, we can cast the eigenvalues into three 
classes, CfjCjjCg. They comprise, respectively, non- 
degenerated, two-fold and more than two-fold degener- 
ated eigenvalues. 

A total of D% = T 1 - 1 + 1 and D% = 2™ eigenvalues 
are found in the classes C™ and C% ■ Eigenvalues in these 
classes do not appear again for any other value of n. 
However, the new eigenvalues remain in the proximity 
of the corresponding ones in the former generation. The 
total of remaining eigenvalues in the class C3 1 is = 
I (3 n_1 — 2" + As the total number of eigenvalues 
increases with 3", the relative number of states in Cf 
and C2 become vanishingly small as n — > 00. 




FIG. 2: Illustration of parity properties for three eigenstates 
at generation n — 2: a) Even state corresponding to the 
largest value E in class Cf; b) and c) Odd states correspond- 
ing to a two- fold degenerate eigenvalue in class Cf. 



Contrary to what is observed with the two first classes, 
eigenvalues in C3 are recurrent, i.e., once they are present 
in the spectrum for a given value of n = q, they will be 
found again for any n > q. Moreover, their degeneracy 
increases with n in a very precise way. Therefore, the 
number of new emerging eigenvalues together with their 
degeneracy can be quantified according to the following 
scheme: the number of new degenerate eigenvalues that 
emerge at generation q is 2 9-3 ; at generation n, the de- 
generacy of any eigenvalue that first appeared at gener- 
ation q is d n>q = - — -q ~ 3 . Note that d q . q — 3 always. 
The first C3 eigenvalues emerging at q = 3 and q — 4 are, 
respectively, and ±\/3- The large degeneracy indicates 
that the characteristic polynomial is highly factorized. 
Nevertheless, this knowledge does not help much in find- 
ing exact roots of the polynomials when q > 5. 



IV. RESULTS 

A. Symmetry properties 

We have scrutinized all eigenstates until n = 9, finding 
a clear relationships among several of their properties 
and the classes to which the corresponding eigenvalue 
belongs. The most direct one refers to the way the states 
reflect the network invariance by 27r/3 rotations around 
the central node. Let us define a parity-like property P 
based on the value of XL =1 a-E,i, so that P — odd (even) 
when this sum is zero (non-zero). We have found that, 
for any P — even state, all sets of three equivalent sites, 
i.e., the three sites that are mapped onto one another by 
rotation of the network by 27r/3 and 47r/3 around the 
central site, share exactly the same amplitude a^. The 
same amplitude is found for all sites that are mapped 
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onto one another by inversion operations over the three 
bisectrices of the original triangle. All C™ states are even 
and, conversely, all CV; and C3 states are odd. 

However, we find two kinds of P = odd states, de- 
pending on whether they belong to the or C3 classes. 
Irrespective of the particular class, the amplitude at the 
central site for any odd state vanishes identically, i.e., 

{a>E,i=N(n-i) — 0)- 

In Figure 2 a-c we illustrate some features of three 
states when n = 2. They correspond, respectively, to 
the largest eigenvalue E7, which always belongs to Cn\ 
and to the second largest two-fold degenerated eigenvalue 
E 5 ,E 6 = l)/2, which is in the class dp. As 

exemplary shown for n = 2, the largest eigenvalue state 
has only positive components (absence of nodes). The 
two states corresponding to E&fi have odd parity and are 
clearly orthogonal. 

The states in are such that, for each set of three 
equivalent sites (by rotation of the network by 27r/3 and 
Att/3 around the central site), the sum of the amplitudes 
a,i is zero. However, no inversion symmetry, as found for 
even states, are observed. 

The states in C3 do not have the same general rotation 
symmetry around the central site, but a large number 
of local 27r/3 rotational symmetries. Indeed, we found 
that, if the corresponding eigenvector appeared for the 
first time in Cf , at generation n it contains N(n — q + 2) 
nodes with zero amplitudes, including the central one. 
For instance, when n = 4, all 12 states corresponding to 
E = 0, which appears for the first time when n = q = 3, 
have iV(3) = 16 sites with vanishing amplitude. The 
3 states corresponding to E = y/3, which appears for 
n = q = 4, have 7 sites with vanishing amplitudes. The 
sites with vanishing amplitudes are exactly those sites 
that are present in the network at generation n — q + 2. 
Further, each set of three sites that appears around these 
zero amplitude sites at generation n — q + 3, have the 
property that the sum to the corresponding amplitudes 
zero. This is illustrated in the Figure 2d, when n = q = 3, 
for a state with eigenvalue E = 0. 



B. Localization properties 

The localization properties of the quantum states were 
evaluated with the help of equation [31 The general fea- 
tures of the results are displayed in Figure 3 where, in 
order to better compare the values of £ for different val- 
ues of n, we draw the value of £(i)/N(n) for individual 
states i. As the spectrum is highly degenerated, the label 
i is chosen to be, in first place, an increasing function of 
the energy. Then, within a degenerated level, the states 
are labeled according to increasing values of £. Since the 
horizontal axis has also been scaled to 1 by dividing i/N, 
the graphs show a recurrent form, with the presence of 
ever fine details as the value of n increases. 

So, despite the shown patterns seem quite irregular, 
it is possible to identify large windows where the value 




FIG. 3: Dependence of £i/N with respect to i/N for n = 
6, 7, 8, and 9. The states are labeled by increasing values of 
E and, within each degenerate level, by increasing values of 






FIG. 4: Absolute value of wave amplitudes for three states 
when n = 7. (a) and (b) correspond to states with E — and, 
respectively, smallest and largest value of £. In (c) we show 
the state with largest value of £, in the class C| and E = \/3. 
The increasing value of £ stays in close correspondence with 
the trend to change from localized to extended character. 



of £/N increases monotonically. These structures cor- 
respond the highly degenerate levels of C3 class, where 
the individual states have been classified according to in- 
creasing values of £. This behavior make it clear that 
the states corresponding to degenerate levels have own 
particular features. 

The behavior of the states corresponding to the classes 
C™ and are much difficult to be perceived in Figure 
3, as they are immersed into the overwhelmingly larger 
number of C3 states. However, the behavior of some of 
them, at specific positions in the spectrum, can be iden- 
tified. For instance, the states corresponding to smallest 
and largest value of E, both in C™ Vn, have very low val- 
ues of £. On the other hand, C™ and states that pre- 
cede the large C3 windows, like those just before E = 
and E = y/3, are characterized by quite large values of 

In Figure 4a and 4b, we plot the absolute value of local 
amplitudes for two n = 7 states: % — 125 and i = 244, 
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ioo. E = 0, largest a - 0.92 



E - 0, smallest a - 0.40 



E = - 3 , largest a - 0.89. 



E =- 3 , smallest a - 0.81 





E - 1.23, a - 1.0 



1=1 , E = smallest eigenvalue, o - 0.05 



i-2, E - second smallest eigenvalue, 



FIG. 6: (a) Dependence on N of the average participation 
ratio (£}, for generations n = 5,6,7 and 8. All states with 
energy E = have been casted into 10 groups, according to 
their values of £. The values of the different slopes, drawn 
in the inset, indicate that most states are extended, (b) The 
participation ratio probability distribution P(^/N) of states 
with E = 0, for n = 7,8, and 9, is largely concentrated on 
large values of £/N for all n. This is in accordance with (a). 



FIG. 5: Participation ratio £ for some specific states, as func- 
tion of N, for generations n = 4, 5, 6 and 7. (a) E = and 
(b) E = — v3 correspond to states of class C3 1 . (c) and (d) 
illustrate localized and extended states belonging classes C™ 
and C2, respectively. 



both with energy E = and, respectively, smallest and 
largest value of £. It is possible to visualize that the state 
with smallest value of £ is characterized by a few large 
amplitude spikes, distinguish themselves from the much 
lower amplitudes in the remaining sites. This picture 
contrasts with that for i = 243, where the amplitudes 
are homogeneously distributed over all sites, evidencing 
properties of an extended state. Figure 4c, for the 
state i = 286, with energy E ~= 1.0495, shows a pattern 
similar to that in 4b, in agreement with the indication of 
an extended state (£/N = 0.4181). 

To characterize the state properties in the n — > 00 
limit, it is necessary to follow the behavior of £ with n. 
As discussed before, truly extended states obey a linear 
dependence between and N(n). However, as the 
number of states increases with the size of the system, 
a criterion must be defined to make the correspondence 
among the states in different generations n. For states 
in C3 , such correspondence is easier to be made, as the 
eigenvalues in one generation will be present in all further 
generations. As this does not happen for eigenvalues in 
C™ and C£ , we will compare them by their relative po- 
sitions to the eigenvalues in the class C3 1 . 

In Figure 5, we show the dependence of £ on the num- 
ber of sites N for some selected states, from which it 
is possible to evaluate the exponent a in £ ~ N a . In 
Figure 5a, the two states with E — and smallest and 
largest values of £ are characterized by a = 0.40 and 0.92, 
respectively. The indication of localized and extended 
character confirm the overall picture shown in Figure 4. 
The result in Figure 5b, for the corresponding states of 
E = — v3, reveal a much smaller variation in the value of 
a £ [0.805, 0.897] for the states corresponding to extreme 
values of £. Similar results are obtained for other values 
of E in C3 , with a clear indication that, despite the oc- 



currence of some localized states, the results in this class 
have a grater tendency of containing states with extended 
character. 

States of the C" class can display both localized and 
extended character as indicated by the sequence of val- 
ues of £ for the largest and smallest eigenvalue and for 
the state with largest value of £. While small slopes 
a = 0.051 obtained for the sequence of smallest eigen- 
values (Figure 6c) indicate localization, the state with 
E ~ 1.239, characterized by a ~ 1, which maximizes the 
value of £ in all generations, shows the opposite extended 
pattern. Similarly, states in class can display both lo- 
calized and extended character, as shown in Figure 5d. In 
the first one, the behavior of £ indicates localized nature 
for the pair of states corresponding to the second small- 
est eigenvalue, while extended properties are obtained for 
the pairs of states with E ~ 0.72. In both cases, small 
vales of a hints at a strong localization character. On the 
other hand, extended properties are found for the series 
of eigenvalues that precede the eigenvalue E = 

Since the relative number of states in classes C" and 
C" decrease with n, let us concentrate on the properties 
of C3 states. A further characterization of their local- 
ization character can be provided by taking the average 
values of £ for the states with the same eigenvalue. For 
this purpose, for each value of n, the states to each de- 
generated eigenvalue have be divided into 10 groups, ac- 
cording increasing values of £. Then, as in Figure 5, we 
draw in Figure 6a (£(n)) as function of N(n). The val- 
ues of a for each subset show, consistently, a tendency 
to accumulate at a large values > 1.0. Another way to 
confirm the result is to draw the the participation ratio 
probability distribution P(^/N) as function of £/N (see 
Figure 6b). For three increasing generations, the posi- 
tion of the peaks of P indicates that the values of £ for 
E = states increase linearly with the size of the system, 
so that their extended nature is preserved in the limit of 
infinite system. These features support our claims of a 
overwhelming extended character of states in the uniform 
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FIG. 7: Dependence of £i/N with respect to i/N for n = 6, 7, 
and 8 for j3 = 1 (a) and /3 = — 1 (b). The states are labeled 
by increasing values of E. In (c), we show the dependence 
of (a) with respect to (3, with a clear indication of a sudden 
change in the nature of the states for /3 = 3. 



model. 



Localization properties of the non-uniform 
model 



Let us now discuss the properties of non uniform model 
mentioned in Section II, where V(i,j) — 1/z". Posi- 
tive values of (3 decrease the relative hopping probability 
from electrons in the largely connected hubs. So, in prin- 
ciple, we should expect a decrease in the electron mo- 
bility and a tendency for localized states. The opposite 
happens for negative values of f3, where the hopping of 
electrons from the hubs is enhanced, strengthening the 
mobility. Since the V(i,j) non-uniformity destroys the 
exact scheme based on three well defined classes, it is 
now better to provide an analysis of the average behav- 
ior of the states, rather than following them individually. 

Our results are summarized in Figure 7. For two dis- 
tinct values of (3 = ±1, we qualitatively see, in Figure 7a 
and 7b, that the dependence of participation ratio £/N 
with respect to i/N has been distorted with respect to 
that of j3 = (Figure 3). Nevertheless, it is still possible 
to identify the position of the (3 = highly degenerated 
C3 eigenvalues. In order to quantify the change in the lo- 
calization properties as function of /3, we draw in Figure 
7c the behavior of the average value of the state slopes 
(a) with respect to (3. The abrupt change in the average 
value of (a) at f3 = 3 hints to an extended-localized tran- 
sition of the wave function at this value. This behavior 
is in accordance with the discussion in the beginning of 
this subsection. 



TABLE I: Dependence of the localization parameter a on dis- 
order strength A. For all investigated values of A, a strong re- 
duction in the value of a in comparison to the uniform model 
is noticed, indicating complete state localization. 



D. Anderson transition 

In this subsection we briefly point out the most im- 
portant features of the eigenstates that result from in- 
troducing disorder into the tight binding model. In 
this case, the on-site energy e; is considered as a ran- 
dom independent variable, describing the local disor- 
der that disturbs the motion of electrons, is consid- 
ered an variable characterized by a probability function 
P(ei) = $(A/2- |ej|)/A, with <I> being the step function. 
The parameter A is a measure for the disorder strength, 
and the main purpose is to identify the possible existence 
of a Anderson transition for some finite value of A. 

As observed in the discussion of the non-uniform 
model, the introduction of the disorder breaks the high 
degree of degeneracy of the C3 states, as well as the 
three classes scheme observed in the energy spectrum of 
the uniform model. Thus, we resort again to a discussion 
of the average properties of the states as a dependency 
of the energy interval. 

The Anderson transition has been studied taking the 
average values of £ for the all states around energy 
\E\ < SE, with SE = 0.001. The results are summa- 
rized in Table 1, where we present obtained values of 
(a) as function of A. The results have been obtained by 
performing M different realizations where, respectively, 
M = 1000, 300 and 50 for n = 6, 7 and 8. 

The results in Table 1 can be compared with those 
obtained for similar models built on complex network 
structures. Indeed, Anderson transition investigations 
on scale-free networks have shown that the localization 
of the electronic states is influenced by the connectivity of 
the network. Thus, the fractal dimension quantifies the 
critical disorder [2CJ, [21| . Our investigation confirms that 
the possibility of the absence of a Anderson transition at 
highly connected networks, as observed in Ref. (2oT |. 



V. CONCLUSIONS 

In this work we presented a very broad investigation of 
the properties of the eigenstates of a tight binding model 
on the Apollonian network. We considered the uniform 
model and two further model versions, where we could 
investigate, respectively, the influence of non uniformity 
induced by the network geometry and by random on-site 
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energies in the eigenstate properties. 

In the regular case, we could establish a precise re- 
lationship between the three different classes in which 
the eigenvalue spectrum can be divided and parity prop- 
erties of the corresponding eigenvectors. We have shown 
that localization property is not related to the eigenvalue 
classes, by the identification of eigenstates, with localized 
and extended properties, within each of the classes. How- 
ever, as the number of states in the class C3 increases 
much faster with n than those of the other two classes, 
we proceeded with a quantification of the average par- 
ticipation ratio of the C3 states. The results indicate 
a scenario in which extended states dominate the global 
behavior of the system. 

By including an explicit dependence of the node de- 
gree on their ability to decrease or increase the electron 
mobility, we have shown that the general state proper- 
ties can change from the quoted extended character into 
a localized one. Our results were based on the ansatz 
that the interaction strength decays with the node de- 



gree according to a power law. It was possible to obtain 
a positive critical value for the change in the property of 
the states. 

Finally, the model has shown to be much more sensitive 
to state localization when random on-site energies are 
assigned to each site. Our results suggest that for any 
non-zero amount of randomness, all states, even those in 
the neighborhood of the band center, assume localized 
nature. 

Physical models constructed on AN's are known to pro- 
duce quite unusual properties. These are due mostly on 
the the existence of large degree nodes as well as with 
the existence of large number of loops. Because of this, 
it is not possible to affirm that other non-uniform disor- 
dered models on AN will not present other non-expected 
features. 
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